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Asymptotic methods are used to describe the nonlinear self-interaction between pairs 
of oblique instability modes that eventually develops when initially linear spatially 
growing instability waves evolve downstream in nominally two-dimensional laminar 
boundary layers. The first nonlinear reaction takes place locally within a so-called 
‘critical layer’, with the flow outside this layer consisting of a locally parallel mean 
flow plus a pair of oblique instability waves which may or may not be accompa- 
nied by an associated plane wave. The amplitudes of these wave's, which are com- 
pletely determined by nonlinear effects within the critical layer, satisfy either a single' 
integro-differential equation or a pair of integro-differential (‘([nations with quadratic 
to quartic-type nonlinearities. The physical implications of these (‘([nations are dis- 
cussed. 


1. Introduction 

Transition to turbulence in boundary layers usually begins with initially linear and 
non-interacting instability wave's that grow to nonlinear amplitudes as they prop- 
agate downstream. The first nonlinear stage of evolution (which might more ap- 
propriately be referred to as a modal-interaction stage) is usually characterized by 
the rapid growth of three-dimensional disturbances due to resonant interactions be- 
tween instability waves and between instability waves and streamwise vortices. This 
phenomenon is usually studied experimentally by artificially exciting the flow with 
small-amplitude noarly-two-dimensional and single-frequency excitation devices. The 
initial unsteady motion just downstream of the excitation device should then have' 
harmonic time dependence and be well described by linear instability theory. In 
most cases, the mean flow is relatively two dimensional and fairly close to a B Iasi us 
profile in the low Mach- number experiments and to the corresponding compressible 
flow in the high Mach-numbcr experiments. (There are actually very few controlled 
excitation experiments at high Mach numbers, but there is some hop(‘ that this 
will change in the near future.) The instability-wave growth rates should therefore 
be small compared to the inverse boundary-layer thickness A~ l at subsonic Mach 
numbers, but can be of the same order as A~ l at sufficiently high supersonic Mach 
numbers (due to the generalized mean-flow inflection point that occurs in this case). 
However, mean-flow divergence effects will usually cause the growth rate to be small 
(relative to A~ y ) bv the time nonlinear effects set in, even in these more unstable 
supersonic flows. This is because', in the latter type of flows, the excitation is usually 
located in the vicinity of the peak local growth rate of the relevant normalized in- 
stability growth-rate versus frequency curve such as the one shown schematically 
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Figure 1. Typical linear growth-rate curve. 


iii figure 1. The growth rate' should therefore decrease as the (constant frequency) 
instability wave propagates downstream into a region where (in most cases) the 
bound ary- layer thickness A will have increased. 

This suggests that the method of matched asymptotic expansions can be used 
to describe these flows: with an inner 1 nonlinear region, in which the instability- 
wave growth rate is small, and a much larger 'outer region in which the unsteady 
flow is governed by linear dynamics, but in which mean-flow divergence effects are 
important (Goldstein k Leib 1988). Once the solutions in these two regions have been 
found, a uniformly valid composite solution that applies everywhere in the linear and 
nonlinear regions can be obtained in one of the usual ways say. by multiplying the 
linear and nonlinear solutions together and then dividing through by their common 
part in the overlap domain (that always exists between the inner and outer regions). 

Smoke wires and other flow-visualization devices can be used to observe 1 the 
flow when the Mach number is sufficiently small. The initial nonlinear (or modal- 
interaction) stage then becomes visible through the appearance of rows of yl-shaped 
structures, which can either be aligned or staggered in alternating rows depending 
on the experimental conditions. The unstaggered arrangement, which was originally 
observed by Klebanoff k Tidstrom (1959) and Klebanoff et al (1962), is usually 
referred to as ‘peak- valley 1 splitting (or K-type transition). It is (as pointed out by 
Kachanov k Levchenko (1984, §5.2)) a complex phenomena that is explainable bv at 
least three different (relatively weak) resonant-type interaction mechanisms each of 
which probably played some role in one or more of the many experiments that have 
been carried out to study this phenomena (see Kachanov et al. 1985; Kachanov 1987: 
Hama & Nutant 1963; Kovasznay et al. 1962; Nishioka et al 1979). A resonant-tvpe 
interaction between oblique instability waves and weak streamwise vortices seems to 
have played a dominant role in the original experiments of Klebanoff k Tidstrom 
(1959) and Klebanoff et al (1962). However, Stewart k Smith (1992) propose a differ- 
ent mechanism that seems to be in good agreement with experimental observations. 

The staggered arrangement, which tends to predominate at the lower excitation 
levels, is usually associated with a weak nonlinearity result ing from a resonant-triad 
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interaction between a pair of oblique subharmonic modes (which, in most cases, orig- 
inated from the background disturbance environment) with the basic fundamental 
two-dimensional mode generated by the excitation device. This type' of interaction 
was originally analysed for Tollmien Schlichting-type (i.e. viscous-type) instabilities 
by Craik (1971) and subsequently by many others who used finite Reynolds-number- 
type approaches and recently by Smith & Stewart (1987) who used asymptotic meth- 
ods. (However, see Khokhlov (1994a) and the discussion section of Wu ct ai (1994).) 
Moreover, there are now a number of carefully controlled experiments (Kachanov 
ct al 1977; Kachanov & Levchenko 1984; Saric & Thomas 1984; Sarie ct al 1984: 
Coke & Mangano 1989) that basically verify the resonant-triad mechanism, but the 
observed growth rates tend to be much smaller than those predicted by the finite 
Reynolds- number theories (Khokhlov 19946). 

Craik (1971) proposed that the resonant-triad interaction could also play a role in 
the aligned or K-type transition, but with the interaction taking place between a pair 
of oblique modes at the forcing frequency and the small two-dimensional instability 
mode that is invariably generated at the first harmonic of this frequency (see ^5.2 
of Kachanov & Levchenko (1984) for a more complete discussion of this issue). All 
of the relevant, modes could then be generated by the excitation device and would 
not have to emanate from the background disturbance environment. However, the 
analysis could not predict the observed gradual transition from a two- to a three- 
dimensional flow structure unless the (common) amplitude of the oblique modes were 
able to exceed that of the two-dimensional fundamental and, consequently, that of the 
(usually much smaller) first harmonic that causes the enhanced growth of the oblique 
modes. This behaviour would obviously be favoured if the oblique modes were' unable 
to suppress the growth of the first harmonic until they themselves became very large 
because this would allow the oblique modes to continue their rapid growth until 
they became larger than the more slowly growing two-dimensional mode generated 
at the forcing frequency. As shown below, the present high-Reynolds-number theory 
actually exhibits this behaviour. In any event, this latter mechanism probably played 
an important role in the recent peak- valley split t ing experiments of Kachanov (1987) 
and Kachanov ct al. (1985) but may not have been very significant in the original 
Klebanoff experiments. 


2. The outer linear flow 

We first consider the initial linear region just downstream of the excitation device 
where tlu 1 instability waves are still small enough so that no significant modal in- 
teractions take place. At supersonic Mach numbers below about 6 or so where 
the so-called first-mode instability is dominant (Mack 1984, 1987). the most rapidly 
growing modes are oblique instability waves, and the first modal interaction to take 
place would probably be the self-interaction between symmetric pairs of oblique in- 
stability waves (Leib fc Lee 1994). In which case, it is appropriate to suppose' that 
the unsteady motion is initiated from a pair of oblique (cqui-amplitude) instability- 
wave modes with the same streamwise wave' number ev r and se*ale'd angular freepiency 
lu'^A/U^ — (\ r c v and equal and e>pposite spanwise wave numbers (±/i). (?/ x is the 1 
characterist.ie: velocity of the flow, and the subscript r is used to denote' the 1 re'al part 
of the' wave' number o anel the phase speed c, as we'll as all othe'r quantity's to which 
it is appended.) These two wave's combine te> form a standing wave' in the' spanwise 
direction that propagates only in the dire'ction of How which is the' situation that 
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most frequently occurs in wave excitation experiments that typically involve longish 
excitation devices placed perpendicular to the How. 

The two-dimensional mode usually t exhibits the most rapid growth at subsonic 
speeds provided, of course, that the mean flow is sufficiently two dimensional. How- 
ever. oven very weak spanwise periodic mean-flow distortions (i.e. streamwise vor- 
tices) can cause appropriate oblique modes to grow faster than the plane wave at 
the high Reynolds numbers being considered herein. This could occur, for example 
(Goldstein k Wundrow 1994). through a kind of resonant-interaction mechanism 
that was first considered by Nay fell (1981) and Nayfeh k Al-Maaitah (1988) and 
later applied to Gdrtler vortices by Hall k Seddougui (1989). In which case it would 
again be appropriate to initiate the unsteady motion from a pair of oblique waves 
of the type described above but with the growth rates equal to the parametric 
growth rates given in Goldstein k Wundrow (1994). But, even when no streamwise 
vortices art' present (or when they are very weak) and the mean flow is effectively 
two dimensional, the oblique inodes can eventually exhibit the most rapid growth 
due to a parametric resonant interaction with the plain 4 wave that exhibits the most, 
rapid growth in the initial linear stage 4 . The oblique modes can then become large 
enough to interact, themselves nonlinearly upon passing t hrough the parametric reso- 
nance interaction stage, which can be treated simultaneously with the self-interaction 
stage if we initiate the unsteady motion from a resonant triad of instability waves 
in the initial linear region a plane 1 fundamental-frequency wave 1 , with scaled angu- 
lar frequency 2uu^A/U. lc . and a pair of oblique equi-amplitude subharmonie waves, 
(again) with the same 1 streamwise wave number and angular frequency, o, and a r c r , 
respectively, and equal but opposite spanwise wave numbers ±/L In this case 1 , the 
term ‘resonance’ implies, among other things, that the three 1 waves all have the same 1 
phase 1 spee'd c r . This occurs (for the small growth rate's and large 1 Reynolds numbers 
that are 1 of interest here) when 

/i = v/3« r , (2.1) 

which me'ans that the oblique instability waves make a 60° angle 1 with the direction 
of flow. We 1 can, e)f course 1 , allow this angle 1 to be 1 arbitrary in flows wIktc an oblique 
mode can grow more 1 rapidly than the plane wave 1 and resonant reaction with the 
latter is not requireel to enhance* the giowth rate of the 1 formeT. 

Our choice e)f the 1 initial linear modes may, at first, se i em somewhat artificial, but 
the 1 line'ar and parametric resonance stages act as narrow band filter's that arc 1 able 1 to 
select out these 1 disturbance's from relatively generic background disturbance fields. 
Me>rea>ver, we eventually show (at the end of §3) that the resonance condition (2.1) 
doe's not have 1 to be satisfied exactly and that the 1 results actually apply to a fairly 
broad range 1 of wave numbers about the resonant condition. 

It is only possible to develop a systematic asymptotic t hea>ry of these phenomena 
when the Reynolds number Re is assumed to be 1 large. Then, since 1 we also require 
that the 1 instability- wave growth rate's be small in the 1 nonlinear region of the flow, 
the initial modal and nonlinear interactions will be confined to a localized region 
eamtreal around the 'critical level’ (Lin 1951), where the mean- flew velocity, sav U c , 
is equal to the common phase 1 velocity c r of the two or three 1 mode's that interact 
the're (sea' figure 2). This occurs because* the 1 flow must be 1 nearly steady in a reference' 
frame' moving downstream with the phase velocity c r at the small growth rate's being 

f Stewart Kr Smith { U)ST) show that non-parallel How effects can sometimes cause 1 the oblique mode's 
to grow faster than the two-dimensional mode'. 
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Figure 2. Asymptotic structure of flow. 

considered herein, and, the corresponding mean flow is then equal to zero at the 
critical layer by definition. There is, therefore, no mean flow to linearize the analysis 
about in this region, and the nonlinear and/or modal interaction effects will then 
come into play at the lowest possible order of approximation there. This explains 
why energy exchange between resonant modes (which share a common critical layer) 
is much more efficient than between lion-resonant mode's. 

The flow outside the ‘critical layer is still governed by linear dynamics, which 
means that it is either (a.) still given by a locally parallel two-dimensional st.reamwise 
mean flow, sav t/ (;</), or perhaps (in the case where the wavelength scale factor a is 
<C 1) a nearly two-dimensional streamwise mean flow, say U(y) Trr 4 t/ 0 (/y, 2Z), plus a 
pair of oblique instability- wave modes or (b) is given by a locally parallel mean flow 
U(y) plus a pair of oblique instability modes accompanied by a first harmonic plane' 
(i.e. two-dimensional) instability wave. The (external) transverse velocity fluctuation 
v is then given (in the general compressible’ case) by 

v = f <m Re [sec 6A(. r„) 4>(y)v' x cos Z + (f /a) 1 /: t4 (l (:r„) «? ( ,(vy)e i2A '] , (2.2) 

wlu're 

X = cm(.r — act ), Z — adz. (2.3) 

c — a c 

ftl = ' r [ f+0 ((TTAv) 

tf<i = (T(f/tr) l/:i x, (2.5) 

the scale factor a ^ 1 has to be inserted in order to simultaneously cover the Q( 1 ) and 
long- wavelength eases, 0 — tan -1 (/3/d), and Re denotes the real part. The stream- 
wise, transverse and spanwise coordinates, normalized by the boundary-layer thick- 
ness zl, are .r, y and z, respectively; t demotes the normalized time, and 9 denotes the 
propagation angle of the oblique mode (which is ecpial to ^tt when the oblique modes 
are resonant with the plane wave, but is otherwise arbitrary). The scaled spanwise 
wave number, streamwise wave number and phase speed /I d and r, respectively, are 
purely real. 

The’ first term on the right-hand side of expiation (2.2) represents the oblique 
modes, while the second term represents the plane wave'. <P and <P {) are the linear 
normal-mode shapes which can, in general, be found bv solving the appropriate 
Rayleigh's expiation (but see below). A and A () , which depend only on the streamwise 
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coordinate (and then only through the scaled streamwise variable x 0 . which varies 
on the length scale of the nonlinear region which, not very surprisingly, turns out to 
be the reciprocal instability- wave growth rate) determine the overall growth of the? 
instability waves and are. therefore, the most important quantities in these equations. 
They are completely determined by the nonlinear dynamics within the critical layer 
and are, in practice, found by equating the velocity jump across the critical layer, as 
calculated from the external linear solution (i.e. the solution to Rayleigh’s equation) 
to the velocity jump calculated from the internal nonlinear solution within the critical 
layer, e and e (e/a) 1 ' 3 are the amplitude scale factors for the oblique and plane waves, 
respectively, where ( is always much less than a. 

Notice that the growth-rate and oblique-mode amplitude scalings o{tjc r) 1/,:S and 
e . , respectively, arc related. This relation ensures that linear growth arid nonlinear 
(or modal interaction) effects will both impact the external linear solution at the 
same asymptotic order. It is dictated by the requirement that the nonlinear stage 
correspond to the first stage of evolution beyond the initial linear region, i.e. that 
the nonlinear solutions match onto the upstream linear solutions in the matched 
asymptotic sense. The Benney Bergeron (1969) parameter A = l/ecr 4 Re, where Re 
is the Reynolds number based on the boundary layer thickness A, is (in the present 
context) a measure of the relative importance of viscous to growth-rate effects within 
the critical layer, i.e. these effects will be of the same order when A = 0( 1). 

The wavelength scale factor a can be set to unity when the initial linear instability 
wave has order-one wavelength as is usually the case in supersonic flow's with 
sufficiently high Mach numbers. Then the linear instability-wave growth rate will be 
0(r ly/:i ) (Goldstein & Choi 1989) as the nonlinear region is approached (which fixes 
the location of this region). 

For the asymptotically more stable flow's, such as subsonic boundary layers with 
sufficiently small ad' verse- pressure gradients (= ()(o 2 )) or with sufficiently w^eak 
spanwise mean-flow' distortions, a wall be small compared to unity, and the lin- 
ear growth rate wall scale like a 4 over most of the unstable region (Goldstein & Lee 
1992: Wundrow et a l 1994; Goldstein & Wundrowr 1994). The nonlinear critical-layer 
effects will therefore come into play over most of the unstable region (and not just 
near the neutral curve) if we take 

<7(f/rr) 1/: * — a 4 , when A = 0(1), cr(e/A cr) 1/M — o\ when A — ► oc. (2.6) 

And for even more stable flow's, such as accelerating boundary layers (Reid 1965; Wu 
1993) with 0(1) pressure gradients, the growth rate will be 0(a 2 ) over the main part 
of the unstable region. In which case, the nonlinear critical- layer effects wall come 
into play in the major portion of this region if we take 

(f/A cr) 1 ^ 4 = a, as A — ► oc. (2.7) 

Finally, w r e note that the phase speeds of the oblique and plane- wave modes, c and 
Co, respectively, will only be equal (i.e. resonance will only occur) if d and 3 satisfy 
(2.1), or equivalently if a and 3 satisfy (2.1) to within order <x(e/( 1 4- A)cr)b /J . 

3. Critical layer dynamics and the amplitude equations 

The lowest-order critical-layer equations turn out to be linear and (in the most gen- 
eral case) correspond to a balance between growth (i.e. non-equilibrium), mean- flow' 
convection and viscous-diffusion effects. The nonlinear and modal interaction effects 
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are weak in the present description, which means that they do not affect the low- 
est order equations, hut enter only through inhomogeneous terms in a higher-order 
problem. This ultimately implies that the scaled oblique-mode amplitude function ,4 
can be determined from a single amplitude equation or that the amplitude functions 
,4 and ,4 0 can be determined from a pair of amplitude equations depending on 
whether or not the parametric resonance stage plays a role in the interaction. In 
the former case, the relevant equation corresponding to the generalized scaling (2.3) 
through (2.5) is given by Goldstein & Choi (1989), Wu et al (1993), Goldstein k 
Wundrow (1994) and Lcib & Lee (1994) as 

= kA + n [ f KsA(x\ )A(x 2 )A*(xi + ;r 2 - x) <l.r 2 cl;/- 1 . (3.1) 

( l' r ./-OC J-o c 

and, in the latter case, are given by Goldstein Lee (1992. 1993), Wu (1992), 
Goldstein (1994) and Mallier L Maslowe (1994) as 


dA(x) 

dx 


— kA(x) + i / K()A 0 (xi )A*(2xi - x) dx { 


J —oo 

*X rX 1 


+ry 


' — OC •/ — OC 


K 1 'IzA(xi)A(x 2 )A* (.7*1 + x 2 - x) dx 2 dxi , (3.2) 


— ^ — KqAq(x) + i fTj f f \K\A[){X\)A{X2) 

dx J-Oc J-OG 

xA*(2x\ + X ‘2 — 2x) + K 2 .4 ( .r i )A 0 (7; 2 ) A * (x i T 2x 2 — 2x)\ dj^daq 

f 

J — OC 

xA*(x\ + x 2 4- 7*3 - 2x) dr 3 dx 2 djq , (3.3) 


' — OC J —o o 


K;iA(xx)A{x 2 )A{x;i) 


where the asterisks denote complex conjugates; x, A and arc suitably renormal- 
ized, and shifted variables corresponding to ;r 0 , A , and A 0 , respectively; and p and 
7 are complex parameters which are dependent on the basic mean flow. The real 
parameter k represents the linear growth rate of the oblique mode, or the resonant 
growth rate of this mode when there is a sufficiently strong spanwise distortion of 
the mean flow. The real part of Kq is the scaled linear growth rate of the plane wave. 
Its imaginary part k,o, k () j, accounts for the initial wave-number shift between the 
oblique and plane- wave modes. 

Notice that these are integro-differential equations of the type first proposed for 
Ross by waves by Hickernell (1984), rather than the usual ordinary differential equa- 
tions of the classical Stuart -Watson Landau (Landau & Lifshitz 1987) theory. The 
integrals arise from upstream history effects that produce a gradual phase shift- 
ing between modes when the nonlinearity takes place within a non-equilibrium (or 
growth-dominated) critical layer. This occurs because the evolution or growth effects 
have a dominant (i.e. first-order) effect on the flow within the critical layer, but only 
weakly affect the flow outside the critical layer. The nonlinear (or wave-interaction) 
terms are therefore influenced by the growth effects when they arc generated within 
the critical layer, but not when they are generated outside the critical layer, as in 
the classical Stuart -Watson Landau theory. 

The nonlinear kernel functions K 0 and K\-K% will be discussed subsequently. 
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They turn out to be simple polynomial functions of the streamwise (and correspond- 
ing integration) variables in the inviscid limit A — + 0, and in the general case involve 
integrals of exponentials and polynomials of the* streamwise coordinates, h'o (explic- 
itly depends on the obliqueness angle 6 and AV/.j = Ku^/n- 

Classical Stuart Watson Landau theory suppresses the critical-layer (effects 
which can only be justified when the Reynolds number is assumed to be sufficiently 
small. For inviscidly unstable boundary layers, this assumption is inconsistent with 
the locally parallel flow approximation (Huerre 1980: Huerrc 1987; Goldstein Leib 
1988; Goldstein jsc Hultgren 1988). For high-Reynolds-number viscously unstable 
boundary layers, classical weakly nonlinear theory is restricted to a rather smallish 
region in the vicinity of the lower branch of the neutral stability curve in which case 
the size of the upstream linear region would have to be extremely small. Moreover, 
nonlinearity usually occurs in the vicinity of the upper branch of the neutral stability 
curve in most of the relevant boundary-layer experiments (Mankbadi et al 1993). 

To be consistent with our requirement that the solutions evolve from an initially 
linear stage, the amplitude equations (3.1) or (3.2) and (3.3) usually have to be 
solved subject to the upstream boundary conditions 

A — a (t V J \ Ao -#• v* ui , as x - -oc. (3.4) 

so that they match onto the linear small growth-rate solution far upstream, or that 
they match onto the appropriate resonant growth-rate solution when spanwise dis- 
tortion effects play a role however, see §0 below for an important exception to this. 
Notice that only the first term on each of the right-hand sides of equations (3.2) and 

(3.3) contributes to these equations when A and Ao are sufficiently small as they 
are initially and that (3.4) is then an exact solution to the resulting equations. 

We include the linear wave-number shift k^ to allow for an appropriate amount of 
wave-number detuning in the analyses, which means that resonance (2.1) does not 
necessarily have to be exact and that the analysis actually applies to a relatively 
broad wave-number range about this resonance condition. And even more generally, 
we could modify equations (3.1) (3.3) to include wave-packet, ('fleets, as in Huerre 
(1980), Smith & Stewart (1987), Smith & Bowles (1992) and Wu et al. (1994) but 
in the interest of simplicity, we do not pursue this issue here. 

When applied, for example, to adverse-pressure-gradient boundary layers, the 
solutions to equations (3.2) and (3.3) are not uniformly valid in frequency as 
A uA /U^_ —> 0. This is because there is a viscous Stokes layer at the wall that even- 
tually contributes a term 

(<TC?U? 

(e/ayripUr^A/U^yri' 

where U[. is the Blasius skin friction, to the scaled linear growth rate k 0 when Auj^ /U^ 
becomes sufficiently small. However, the relevant solutions can easily be made uni- 
formly valid for all frequencies (except in the immediat e 1 vicinity of the lower branch) 
by simply replacing the relevant linear growth rate's (k and k () in equations (3.2) and 

(3.3) . respectively) by 


- - . 4 [cufUl ' 

(CTC)>U? 


(3.5) 
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4. The mean-flow change 

A significant feature* of t lit' present high-Reynolels-numbe*r approach is tliat t lit* 
nonline'ar critical-layer interaction produce's a spa mvise- variable mean-flow change' 

ii = ( Mi' fi u { !)..v {) )v 1 ' z . (4.1) 

tliat is of the same order as the oblique-mode instability wave' (see' expiation (2.2)) 
that initially produces the 1 interact ion. However, the* associated cross- flow velocities 

r — fT(((/er) ]/ i Re* c ( ,('“ >lZ . tr = f{t/rr) l/ ' Re a\\e lxZ (4.2) 

turn out to be somewhat smaller than this. 

In the remaineler of the paper, we* discuss the* implications of the* fundamental 
equat ions (3. 1 ) (3.3). 

5. The pure oblique-mode interaction 

We begin by conside'ring the* case* where only the* oblique mode's ('liter into the* 
int eu act ion. We have* se*e*n that this situation is redevant to su])(*rsonic boundary 
layers and to subsonic flows with sufficie'iit ly strong streamwise' vort ice's. It should 
therefore* be* highly r<*le*vant to the* original Klebanoff et ai (1902) experiment- (i.e. 
K - 1 y ] >e 1 t rausil ion). 

T he* ('onnnon eiblicjue'-mode* amplitude* is now ele'terminexl by exjuatiem (3.1). Its 
keriie’l functiem l\o is redativedy ce>m])licat('d whe*n viscous ('fleets are* retainexl. as in 
Wu ef al (1993). but in the* invisciel limit originally (‘onsiderexl by Golelstein A r Choi 
( 1989) and Gedelstein A' VVundrow (1994). it is simply 

f<„ - y an" 0 cos 20(.r — .r \ )[(./' — ,r\) 2 + (.r — .r-j) 2 — cos 2ff(.v — .('j )(.(', -.r 2 )], (5.1) 

provided that rr 1 in the* e’ase* where* strexunwise vort ice's play a role. Wunelrow 
A r Goldstedu (1994) also consideTcxl the* stre*amwise* vortex problem in the e)rde*r-on(* 
wa.ve*lengt h. i.e*. t he* n = 1 limit. The'v show t hat t he* instability-wave* amplit ude' is 
still given by (3.1), providexi the* nemline'ar kernel functions is take*n to be a slight 
generalization of (5. 1 ). 

I\(t vanishe's whe*n 0 — anel the invisciel solution to (3.1) de»velops a singularity 
at a finite* denvnst re*am position (Goldste'in Choi 1989), say ./* s , at all othe'r angle's. 
.4 the'refon* e'xhibits e'xplosive* growth at .r s . wit h t he* local asymptotic behavior be*ing 
give'll by (Goldste'in Ohe>i 1989; Shukhman 1991) 

ooy- -*-*• (r, - 2) 

whe*re* the* re*al parame'te'rs a and (j) are' re'latexl to the* original param<*t-e*rs k and y 
through (jiuulrat lire's. Figure* 3 is a plot of the* scalexl amplitude' function veTsus the 
sealed stre'amwise' coordinate', as calculatexl nume'rically from ('([nations (3.1) and 
(3.1) for n 1.2 and various value’s of 0. The* curve's show that the solution initially 
follows the* linear growth corresponding to 


and that the* e'xplosive* growth oceans very suddenly once nonlinearity come's into 
play. The * dashed curve's are* the* local asymptotic expansions calculated from expia- 
tion (3.2). This re'sult implies that the* ove*rall wave'-nuiiibe*r/growt li-rate scaling is 
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Figure 3. A plot of the sealed amplitude against the scaled streamwise coordinate: (a) 0 — 15°; 
(b) 0 ~ 30°: (e) 0 - (j(F: (<l) 0 — 75°. Solid lint*s: numerical solutions; dotted lines: local 
asymptotic solutions (from Wu cf al. 1993). 

preserved right up to the singularity when a — 1 , which means that the* overall 
asymptotic structure remains intact until the instability-wave amplitude becomes 
0(1) everywhere in the How. and that the motion is then governed by the full non- 
linear Euler equations in the next stage of (‘volution. 

However, the growth-rate amplitude scaling is not preserved in the long wave- 
length limit (7 — > () (corresponding to. say, the weak streamwise vortex-amplification 
mechanism). In this ease, the critical layer expands to fill the inviseid wall layer that 
surrounds the critical layer, causing the How to become* fully nonlinear while the* 
instability amplitudes are still small. The next stage of (‘volution is then character- 
ized by a three-layer structure and is governed by the three-dimensional unsteady 
'triple-deck' equations, but without the viscous terms (Goldstein k Lee 1992). This 
does not, however, imply that the relevant scaling is the usual triple-deck scaling in 
this stage. 

Wu at al (1993) showed that explosive growth also occurs in the viscous case 
and that the local asymptotic behaviour in the vicinity of tin* singularity is still 
given by (5.2). However, they also showed that (as in Goldstein k Leib (1989) and 
Leib (1991) ) then* is a certain range 1 of parameters where explosive growth does not 
occur when the viscous parameter A exceeds a certain (usually very large) value. The 
instability wave will then reach a peak amplitude at souk* fixed streamwise location 
and subsequently undergo viscous decay downstream of that point. 

6. The parametric resonance interaction 

Now suppose that the scaled oblique-mode amplitude /I is very small and remains 
that way during the entire resonant, interaction. Notice that this includes the case 
.4 — 0 (f/( 7 ) 1//: \ where the oblique mode has the same amplitude scaling as the plane 
wave (as was originally pointed out by Goldstein k Lee 1992). 

The last term can be neglected on tin* right-hand side of equation (3.2). which 
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then heroines 


— — — k.A + i 
dx 


K() A Q ( x i ) .4 * ( 2 x i - .r) • 


( 6 . 1 ) 


while? the plane- wave amplitude equation (3.3) reduces to the linear growth-rate 
equation 


di () _ _ - 

— — — K()/i(), 
f r 


(6.2) 


which merely reflects the fact that there is no hack reaction of the oblique mode 
on the plane wave. It may seem rather surprising that this occurs even when the 
oblique-mode amplitude is much larger than that of the plane wave, but the critical- 
layer velocity jump that would produce back reaction at this level turns out to be 
identically zero. It is worth noting that the back-reaction effects would have to be 
quadratic in the oblique-mode amplitudes if they occurred at the equi-amplitude 


stage. 

Since the second member of the oblique-mode equation (6.1) is now linear in .4. 
we refer to it as the parametric resonance term. Its kernel function, which was first 
given in Goldstein fc Lee (1993), is 


K„ = §(x - 2 a) 2 exp(-§A(.r - ) 3 ), 


(6.3) 


where A is a suitably renormalized parameter corresponding to A. Goldstein & Lee 
(1992) give an analytical solution to (6.1) through (6.3) for the inviscid limit A — 0 
and Wundrow ct. al (1993) extend it to the viscous case where A = 0(1). These 
solutions show that the oblique-mode instability-wave amplitude can be represented 
by an infinite series of terms each of which exhibits exponential growth. They also 
show that .4 tends to be dominated by the lower-order terms at small values of i\ 
but that the higher modes rapidly come into play and the infinite tail' of the series 
eventually determines the behaviour of the solution at large values of x. This leads 
to the conclusion that 


(■() exp( ±i arg L4 0 ) exp(K 0r f/5) exp 


(ii 0 ) 1/4 di 


as x 


oc. 


( 6 . 4 ) 


provided that the shifting of the coordinate is correct to 0(a) in the long-wavelength 
limit, where a « 1 and k — | k 0 . Here, ,r 0 is a shifted coordinate corresponding to 
j\ c { ) is a real constant, and .4 ( , is given by equation (6.2). 

Notice that K {) (as given by equation (6.3)) becomes highly concentrated around 
x = X\ in the strongly viscous limit: 

A — > oc, with k = A 1//3 k = 0(1). (6.5) 


Equation (6.1) therefore reduces to the ordinary differential equation 

^ = kA + (6.6) 

dx 

where 

x — A _l//: *x, A — A _1/M .4, and A 0 (x) = A _2/3 A () . (6.7) 

The limit (6.5) corresponds to (among other things) the flat plate or Blasius boundary 
layer, i.e. the flow in which the resonant-triad interaction was first analysed by Craik 
(1971). In fact, equation (6.6) is within a constant factor of the equation obtained by 
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Graik (1971). wlio used convent ional Stuart Watson Landau lht*ory (Stuart I960: 
Watson 1 !)(>(): Landau Mishit/ 1987) tugcth<*r wil h finite Rewnolels-numbeT-typc 
arguments to <le*rive his result. Th<' corresponding limiting form of the general plane- 
wave amplitude (‘(illation (AM) is still the linear equation ((72). 

Virtually all suhhannonie transition experiments have been carried out in flat- 
plate boundary layers with very- low frec-stream Mach numbers so that the Blasius 
boundary-layer solution provides an appropriate description of the mean flow. The 
ainplitiide/growlh-rate sealing for linear instability waves in the major portion of the 
unstable Reynolds-number range is then given by equation (2 .(i) and as shown, for 
example*. In Bodonyi A: Smith (1981) and Goldstein cf. al. (1986). R< — 0(fx“ m ) and 
2a; T A/G v = (){fT 2 ) in this range*. In whie'h case, it follows from equation (2.6) and the 
definition of A that e -- err 1 " and A — rr~' i/2 t and. therefore (in vi(‘w of (‘filiation (0.5)), 
that both terms on the right-hand side of (‘filiation (3.5) are of the same order. This 
scaling is valid at all sufficiently large Reynolds numbers including those correspond- 
ing to the upper branch of the neutral stability curve. Mankbadi cf al (1993) point 
out that the initial parametric resonant interaction first becomes significant in the 
vicinity of (and in souk* cases downstream of) t he* upper branch in virtually all sub- 
harmonic transition experiments carried out to date. They also note that it occurs at 
relatively small values of the frequency parameter 2W A l \ R( = a 12 presumably 
because, as theory suggests, the relative strength of this interaction increases with 
decreasing a. The present asymptotic theory, which holds at small values of dr. should 
therefore provide 4 a reasonable description of this phenomenon. 

The classical high-Reyuolds-number asymptotic solution for the upper branch of 
the neutral stability curve (Reid 1966) is the relevant solution in the initial linear 
and non-interacting si age. This solution is basically inviscid except, for a thin viscous 
wall layer and a relatively thin critical layer which is asymptotically distinct from the 
wall layer. The lowest -order critical-layer equation corresponds to a balance 4 between 
mean-flow conven t ion and viscous diffusion effects. 

As in the general case discussed above, the first modal interactions occur within 
the critical layer once the linear plane- wave amplitude becomes sufficiently large', but 
since* the lowest-older critical-layer equations now correspond to a balance* between) 
linear conve'ction and viscous diffusion e*flex“ts, tlu* oblique* instability-wave amplit ude' 
is dctcrminexl by the* viscous critical-layen* amplituele equation (6.6). However, this 
equal iem. together with equatiem (6.2), shows that the* oblique-mexle growth rate* 
(d.l/el.r)//l continue's to increase as the instability waves propagate elowmstream. so 
that the* non-equilibrium (e>r growth) effects. which are' missing in the* lenvest-ordcr 
crit ie*al-lay(*r expiations, must e*ve*ntually come* into play causing expiation (6.6) le) 
bexxmie* invalid. 

In fact. expiatiems (6.2). (6.6) anel (6.7) imply that ((Taik 1971: Wundrow ci al. 
1993) 

A ~ Go exp(i7r/ 1 ) e*xj) ( k.r H vxp(ko r .r) ] . as ,'r — * x. (6.8) 

V 4K.« r J 

where* ( f () is a n*al constant: we* have* clmsexi tlu* origin of the <x)ordinat.es so that 

A {1 = exp (£<),■./*), (6.9) 

anel, for simplicity, we* assume* that k is re*al. 

Notice that expiation (6.8) does ne>t reduex* te> the limiting form e>f expiation (6.-1) as 
A — ► x,. which me*ans that the* limits A — > x and .r — » x cannot be* intcrevhangexl anel. 
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consequently, that there must be souk 1 intermediate solution that connects asymp- 
totic solutions (G.4) and (G.8). In fact, Wundrow et al (1993) show that approxima- 
tion (6.6) becomes invalid when kx = 0(ln A 2 ' 3 ) and that the non-equilibrium ('fleets 
become of the same order as the viscous effects for larger values of x . at which point 
the How begins to evolve on the faster scale 1 

.7 = A 1/:1 (4 - In = A 2/:5 .7 i , ((>. 10) 

where 

x, =x- T— lnA 2/: * (0.1 1) 

K() r 

is an appropriately shifted coordinate on the ./‘-scale, and the oblique-mode amplitude 
is determined by the fully non-equilibrium equation (G.l), but with /lo(i') treated as 
a slowly varying function of x and the linear growth term treated as a higher-order 
effect. The relevant solution has the yvkb.j form (Goldstein 1994) 

.4 = Co exp(iir/4) oxp(A 2/:! ) ‘ b(£) d£. (6. 12) 

where t he prime denotes differentiation with respect to ;/\ q — k/k 0 , b is determined 
by the transcendental equation 

ifioOr,) / cxp(-2c , -26(i,)C)< a aC. (6.13) 

Jo 


A = -A, 
A 


and C> is a real constant. 


(6.14) 


A G (J (^)V^nr exp(i7r/4) exp + - ~z — (exp(K 0r ,/T ) - 1 )j , as x { -► -oc, 

(6.15) 

which nu'ans that the solution (G.l 2) will match onto the asymptotic expansion (6.8) 
of the solution to the viscous-limit equation (6.6) if we take 


Notice that Ao — ► 0, b —> and, consequently, that 

3 A 2 /3 


C) = \Z«0rCl(|A 1/(i ) f, CXp ( —^4 j • 

v 4 Ko, / 


(6.16) 


This shows that the critical-layer dynamics are eventually controlled by non- 
equilibrium (or growth-rate) effects, even in the Blasius boundary layer, and that 
the uniformly valid solution for the instability-wave amplitude is ultimately deter- 
mined bv non-equilibrium equation (6.1) and not by the viscous-limit equation (6.G). 

Figure 4 is a plot of the oblique-mode growth rate as calculated from the full 
non-equilibrium equation (6.1) for various values of the nearly constant scaled plane- 
wave amplitude A {) . The straight line is the result obtained from viscous limit equa- 
tion (G.G). Notice that the non-equilibrium growth rates are consistently lower then 
the corresponding equilibrium values, which might, as pointed out by Khokhlov 
(1994), just explain the discrepancy between the 1 growth rates predicted by the clas- 
sical (Talk-type theories and those observed experimentally. 
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Figure 1. Oblique mode growt h rate as a function of scaled plane-wave amplitude. (Prepared 

by Dr Sang Soo Lee.) 

7. The full resonant-triad interaction 

Equations (6.12) and (C.13) show that the oblique mode continues to grow (when 
A{\ is given bv equation (6.2)) and must therefore eventually become large enough to 
not only react back on the plane wave but also interact nonlinearly with itself. The 
plane wave and oblique mode will then both evolve on the faster scale x as defined 
in equation (6.10). 

The simplest way to show this is to notice that the viscous parameter A can 
l>e scaled out of the general equations (3.2) and (3.3) by introducing the scaled 
dependent and independent variables (6.10), (6.14), along with Aq = Ay j A^, and 
then replacing the* linear growth rates k and k 0 by the scaled growth rates kj A 1/3 and 
K () /A 1//; b respectively. Then, aside from the vanishing of the linear growth terms, the 
resulting equations will remain unchanged in the limit A — > oc, with k (as defined 
by equation (6.5)) and the barred variables held fixed. These latter equations do 
not possess solutions that satisfy the linear upstream boundary conditions (3.4), but 
they do possess solutions that satisfy the alternative conditions 

.4 — * « (0) exp(i7r/4) exp(6 0 x), ,4 0 — > 1, as x — — oc, (7.1) 

where 

bi>='if C 2 exp(-C' , -2fe„C)dC, (7.2) 

Jo 

and therefore match onto the x\ — + 0 limit of the parametric resonance solution (6.12) 
and (6.13) and the linear plane- wave solution Aq = exp(k {) xi). 

The previous results show that these latter solutions match onto an intermediate 
viscous parametric resonance stage (which is governed by equation (6.6)) and, con- 
sequently, onto the same upstream boundary conditions as equations (6.9) and (6.6) 
(i.e. equation (3.4) with the ~ replaced by ") provided, of course, that 

a m = 0(A" /6 exp(-3A 2/3 /4Kor))- 

This means that the fully interactive stage will be governed by the full non- 
equilibrium equations (3.2) and (3.3) if the oblique modes are exponentially smaller 
than the plane wave at the start of resonance even in the Blasius boundary layer. 
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The smaller linear growth rate of the oblique modes could easily cause this situation 
to occur even when all the modes have the same amplitude at the start, of the 
linear stage. 

We now return to the general resonant- triad equations (3.2) and (3.3). The signifi- 
cance of the various terms in equation (3.2) has already been discussed. The relevant 
kernel functions are given by (5.1) and (6.3). However, we have not, as yet, discussed 
the nonlinear terms in equation (3.3). They account for the back reaction of the 
oblique mode on the plane wave with the first group representing a kind of mutual 
interaction. The relevant kernel functions are given by Goldstein & Lee (1992) and 
Mallier & Maslowe (1994) in the inviscid limit, and by Wu (1994) in the general case. 

The last term in equation (3.3), which is quartie in the oblique-mode amplitudes, 
does not involve the plane-wave amplitude at all. The early, i.e. finite Reynolds- 
number-type, analyses of the resonant-triad interaction (see, for example, CTaik 
(1971) and other references) involve a corresponding back-reaction term that is. 
however, only quadratic in the oblique mode amplitudes. The kernel function for 
this last term (of equation (3.3)) is also given by Goldstein & Lee (1992) and Mallier 
& Maslowe (1994) in the inviscid limit and bv Wu (1994) in the general case. 

8. Concluding remarks 

In most boundary layer flows, it is the oblique-mode instability waves that ul- 
timately exhibit the most rapid growth either directly from the initial linear 
stage or indirectly through an intermediate parametric resonance stage. The cu- 
bic self-interaction between the oblique-mode instability waves is one of the first 
strictly nonlinear interactions to come into play as the instability waves evolve down- 
stream in such flows. This interaction can have a dominant effect on the subsequent 
instability-wave development producing a local singularity (and consequently ex- 
plosive growth) at a finite downstream position in the inviscid limit and sometimes 
producing viscous decay when viscosity is present (Goldstein 1994: Wu 1993). 

The more or less general case is described by equation (3.1), or by equations (3.2) 
and (3.3). but depending on the initial amplitude ratio and the external parameters, 
various limiting forms of these equations can apply to different regions of t he flow 
giving rise to a wide variety of different phenomena. The nonlinear interaction also 
produces a spanwise variable mean-flow change in the linear region outside the critical 
layer. It is of the same order as the oblique-mode instability waves in the inviscid 
case, but can be even larger than these in the strongly viscous case' (Goldstein 1994). 

The author thanks Professor R. E. Kelly for first, suggesting that three-dimensional nonlinear 
critical layers would be scientifically interesting, Professor F. T. Smith for suggesting that the 
resonant-triad analysis be continued into the fully coupled stage, and his colleague. Dr Sang Soo 
Lee, for preparing figure 4. The author also thanks his colleagues, Dr Sang Soo Lee, Dr David 
Wundrow, Dr Lennart Hultgren, Dr Reda Mankbadi and Dr Stewart Leib for their helpful 
comments during the course of this work. 
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